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ABSTRACT 


A method is developed for computing the ground temperature 
accurately over both the diurnal and annual cycles. The ground 
is divided vertically into only two or three slabs, resulting in 
very efficient computation. Seasonal storage and release of 
heat is incorporated, and thus the method is well suited for 
use in climate models. 

Intvoduation. The ground temperature should be computed a^c- 
curately over both the diurnal cycle and annual cycle in a cliimte model. 
The diurnal and seasonal cycles are important because of the highly non 
linear way in which latent and sensible heat fluxes from the surface de- 
pend upon temperature. Seasonal storage and release of heat is aiso Im- 
portant because it can amount to several watts per square meter and can 
si gn ificantly influence the ground temperature (Taylor, 1976). 

The basic problem is to solve the one-dimensional heat conduc- 

tiou equation 
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where c is the heat capacity of the ground per unit volume, X is the 
thermal^ conductivity and K is the thermal diffusivity , subject to num- 
erically specified heat flux ^t the upper boundary; 


= F - F, - F, - F . (2) 
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F is the absorbed solar radiation, F, the net thermal radiation, and 
F are the sensible and latent h^x fluxes from the ground to the 

afeiosphel^. , ... . 

A straightforward finite difference approach with computation oi 
temperatures at several subsurface depths can yield accurate results. 
However, that method has the disadvantage of requiring (1) computation 
and storage of temperatures at several depths, and (2) a short time step, 
of the order of a few minutes or less. 

Slab model. An approach requiring less resources, is to divide 
the ground into a small number of slabs and solve for the mean tempera- 
ture of each slab. In this method we use the average temperature of the 
upper slab to simulate the surface temperature. That layer must be suf- 
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diurnal temperature response Is representative of 
that for the upper skin of the ground; on the other hand, it should be 
thick enough to allow long time steps. The minimum number of ground 
layers that are needed for a climate model is two, since a deep lower 
layer is needed for seasonal heat storage. 

*« *1. Jlif our slab model are derived under the assump- 

* capacity and conductivity are uniform in each slab, and 

wira the temperature profile within each slab approximated as a quad- 
ratic function of depth: 

Tj(z) = a ^ +8^z + Yj. (3) 


For ^o layers we obtain the alx parameters from the six equations: 
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^ = -h~^2 l«iterface. The heat 
crossing the z = -z^ Interface Is: ^ ^ 
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The final equations for the mean temperatures of the two slabs are: 

df. 
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“dt = ^(0) - F(-z^) 
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dTg 

V2“dt = • (7) 

We use a simQar formulation for ocean ice but the no flux condi- 
tlon at z = "Z^-Zg is replaced by: 

T2<-z^- z^) = (8) 

7^5k***® freezl^ point for ocean, is taken as 271. 6 K. The heat cross- 
ing the z = interface is 

F(-. , . ‘^l-18f,.6T^-2F(0)z,A, 

4Zj/Aj + aZj/Xj 
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and the heat crossing the ice-ocean interface: 


F(-Zi - z„) = 




2 


The final equations for the two layer model for ocean ice are: 
dT 

^i°i“dt = • ai) 

dT2 

V2-dT " F(-z^)-F(-Zj- V . (L2) 


The equations for our three layer model are determined in a similar 
way. 

The thicknesses of the layers were determined empirically by 
comparison with an exact periodic solution and an accurate 21 layer fin- 
ite difference solution. For Kg = 4 x lO'^nr s“l tiie optimum depths for 
two layers were found to be z,<lo cm, z 'v. 4 m, while for three layers 
5 cm, Zg^ 25 cm, z„'V4m, For other values of Kg these depths 
should be scaled in proportion to Ihe square root of Note that for 
the particular thermal dlfifusivity which we employ th# damping depth, at 
which the amplitode of the temperature variation is reduced by e"l(cf. 
Sellers, 1972), is''' 10 cm for the diurnal temperature wave and '^2 m for 
the annual temperature wave. 

To illustrate the errors in the two and three layer models , we 
show the results of computations for the case Kg = 4 x 10“^2g-l 
same as case 1 of Deardorff, 1978). For these fests we considered the 
simple case of no atmosphere above the ground (Fj^=Fg=0), with 


^sw ^"^g^ 
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gromd albedo A = 0. 24, solar constant 1. 354 x 10^ W and ground 
emissivity e = 0.^; 0 is the solar zenith angle and a the Stefan-Boltz- 
mann constant. This idealized case yields unrealistic temperature vari- 
ations at tile surface, but it simplifies the computations and their inter- 
pretation. 

Fig. 1 compares the diurnal and annual surfece temperature vari- 
ations computed with several different models. The computations are 
for the simplified case of no atmosphere described above. Our two layer 
and three layer results are about equally accurate for the annual temper- 
ature variation, as eiqiected since both models use one deep layer for an- 
nual heat storage. Three layers are clearly needed if it is desired to re- 
present both the amplitude and phase of the diurnal ground temjierature 
variations, as well as the seasonal heat storage. We emphasize that the 
Airakawa (1972) and Deardorff(1978) schemes were not designed for use in 
climate models, and we have included them only to illustrate the necess- 
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Fig.l Diurnal and annual surface temperature variations computed 
with several different models for the simple case of no atmosphere 
and thermal dlffusivity Kg = 4 x 10’'7m2s‘i. The Arakawa (1972) 
and Deardorff (1978) schemes were not designed for use in climate 
models, and are included here only to Illustrate the necessity of 
having a deep layer for annual heat storage. 

ily of having a deep layer for annual heat storage. Arakawa scheme 
has one layer which is used to represent the diurnal temperature varia- 
tion. Deardorff’s method has two layers, with the lower layer damping 
the upper layer so as to give a realistic phase to the diurnal temperature 
variation. 

Fig. 2 compares the diumally averaged heat flux into the ground 
for the same models as in Fig. 1. The Arakawa (1972) model has a ne- 
gligible diurnally averaged flux, since it consists of a single thin layer 
with a zero flux lower boundary condition. Although the two layers in the 
Deardorff (1978) model are thin, there can be a substantial diumally ave- 
raged heat flux into (or out of) the ground because the method does not 
employ a zero flux lower boundary condition. 

Di. sous si on. The ground temperature model we have presented 
permits efficient and accurate computation of both diurnal and seasonal 
temperature variations with a minimum number of layers. It is there- 
fore well suited for use in global climate models. 

The diumally averaged heat flux into the ground is usxially small 
in comparison to typical radiative, latent and sensible heat fluxes, even 
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r ig. 2 Diumaiiy averaged heat flux into the ground for the same 
models as in Fig. 1. The cjalculations are for the simple case of 
no atmosphere at 45 N latitude, as described in the text. Arakawa’s 
model, since it is a single lliin layer with zero flux lower boundary 
condition, has negligible diumaiiy averaged heat flux. 

for this idealized case in which the moderating effects of the atmosphere 
are excluded. Thus a large impact of seasonal heat storage on short 
time scales should not be expected, but the cumulative effect on season- 
al time scales may be significant. Specific tests of the impact will be 
included in the tiiree-dimensional climate model described elsewhere in 
these proceedings. ; A previous test has been made in a general circu- 
lation model by Delsol et al. (1971), who used a ground thickness of 5 m 
with a lower boundary condition of T=280 K; however, they ran the model 
only 14 days which is too short for an adequate test. 
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